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Error Analysis of Phase-Integral Methods. 
II Application to Wave-Penetration Problems 

F. W. J. Olver 

(July 1, 1965) 

A study is made of the differential equation 

d 2 wldz 2 ={j{z) + h(z)}w, 

in which /(z) and h{z) are real and regular on the real axis, and/U) has exactly two zeros there. Strict 
error bounds are derived for the coefficients in the formulas which connect the asymptotic solutions 
at z = + o° with the asymptotic solutions atz = — ». Applications are made to two physical problems. 

Key Words: Wave penetration, connection formulas, error bounds, potential barrier, transmis- 
sion coefficient, harmonic oscillator. 

1 . Introduction and Summary 

The purpose of this paper is to illustrate the theory of the preeeding paper [l] 1 with two 
problems arising from the one-dimensional wave equation 

— ={Az) + h(z)}w. (1.01) 

In both problems the functions f(z) and h(z) are real when z is real, and in the region considered 
f(z) has no singularities and just two zeros, both of which are real and simple. The function 
h(z) is to be regarded in some sense as being small compared with/(z), but we do not need to formu- 
late this assumption precisely at the outset. 

In the first problem (sees. 2, 3)/(z) is positive between the zeros. The object is to determine 
the ratio of the amplitude of an oscillatory solution at z = -h°° to the amplitudes of its components 
at z = — oo. 

In the second problem (sees. 4, 5) f(z) is negative between the zeros. The object here is to 
determine eigenconditions admitting a solution having subdominant (exponentially small) character 
at both — oo and + oo. 

Both problems have been treated frequently in the literature; see the references cited in [1], 
and [2] — [7] given at the end of this paper. Generally, the viewpoint adopted has been that /(z) 
contains a multiplicative factor u 2 , h(z) is independent of u, and approximate solutions are sought 
which are asymptotically correct for large values of the parameter u. In the present paper, we 
establish explicit strict error bounds for the approximate solutions. Asymptotic properties with 
respect to parameters are immediate consequences of the forms of the bounds. 

Two methods for solution were given in [1], We shall concentrate on the first (Zwaan's 
method). As we shall see, from the standpoint of deriving error bounds this affords the simpler 
approach to the present problems. 

2. The Overdense Potential Barrier 

In this section we suppose that f(z) and h(z) are regular in an unbounded, simply connected 
complex domain D, which includes the whole of the real axis, and that the only zeros of f(z) in 
D are simple zeros at z~±a (a> 0). When —a<z< a,f(z) is positive. 

1 Figures in brackets indicate the literature references at the end of this paper. 
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Let us consider the turning point at z = — a. In order to apply the theory of [1], section 5, 
we have to exclude the other zero z= a from the domain of consideration, and we do this by intro- 
ducing a cut sd along the real axis from z=a to z = -h°°. We shall refer to the upper and lower 
sides of this cut as s^i and j/ 2 , respectively. As in [1], section 5, we define 

f(z)= r {f{i)Y t2 dt. (2.01) 

J -a 

The associated principal curves^- and subdomains D, are indicated in figure 1. *% \ is the part 
of the real axis which extends from — a to — ». ^ 2 ,^ 3 are conjugate curves emerging from — a at 
directions + 7r/3. They cannot intersect the real axis at any other point. We define £/(z) to be 
the branch of £(z) which is continuous in D cut along s$ and ^j, and has positive real part in D,. 
Similarly, g\{z), #2(2), gz(z) are multiples of any branches of g(z) = {f(z)}~ 114 which are continuous 
in D cut along s& and the corresponding^, and satisfy eq (5.08) of [1]. 

As in [1], let q be points at infinity on arbitrary ^-progressive curves J^j lying within D, such 
that |Re f(c/)| =°o. £?i is taken in the upper half-plane, and if 2 , if 3 are complex conjugates, but 
no boundary of D may intervene between if 2 and if 3 in the neighborhood of infinity. The other 
assumptions we make are 

(i) The variation of the function 



F = 



j\^-^S(^)} dz (2 - o2) 



converges as z^> Cj along ifj, and also as z^±°° along the real axis. 

(ii) ^-progressive curves can be found in D linking: (a) Ci with c 2 ; (b) C\ with c 3 ; (c) c 2 or c 3 
with any point of the interval (— °°, — a); (d) C\ with any point of sd\. 

(hi) 6 2 (c 3 ) = 0; e 2 (z) and e 3 (z)->0 as z— > — °°; €1(2) — » as z^ + 00 onij. 

(A brief discussion of the determination of ^-progressive paths is given in the remark at the 
end of [1], sec. 5, and sufficient conditions for the fulfillment of (hi) are given in sec. 7 of the same 
paper.) 

With these conditions, the fundamental connection formula ([1], sec. 5) for the solutions 

Wj{z) = gj(z)e-¥ Z) {1 + €j(z)} ( 2.03) 

of the differential eq (1.01), becomes 

wi(z) = e^ 3 {l + ei(c 3 )}w 2 (z) + e~^{ 1 + €i(e 2 ) W*). (2.04) 
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We fix the Wj{z) by specifying that g\(z) is the branch of {f{z)}~ 114 which takes its principal value 
on the join of — a and +a, and interpret the branches of the two sides of (2.04) in real form, as 
follows. 

When z E .stf h we have 



2.05) 



fi(z) = a-i f ^{fV' 2 dt, gl ( z ) = ^/4|/|-i/4 5 

J a 

where 

«=|"j/|" 2 ^ (2.06) 

and is positive. Therefore from (2.03) and hypothesis (iii) 

M ; 1 (z) = e-«|/|- 1 / 4 {l + o(l)} exp(; (' \f\ ll2 dt + ^ni) (z^ + oc). (2 .()7) 

This expression represents a transmitted wave ([8], p. 1100; [6]). 
When — °° < z < — a, we have 



6« = i f ° \f\ ll2 dt, g2(z)=e-™lu\f\- 



1/4 



(2.08) 



H 



ence 



e 7 "/ :i {l + e,(c,)}^2U)=|/|- 1 / 4 {l + € 1 (c.3) + o(l)} exp(-i [ " \f\^dt+^m 



(2.09) 



This is the incident wave, moving in the same direction as the transmitted wave. Similarly 

1 

4 



e-**/3{l + €i(c 2 )}w 3 (2:)= |/|- ,/4 {l + €i(c«) + o(l)} exp^ i J ° |/| 1/2 A— jwi) 



(2.10) 

This is the reflected wave, moving in the opposite direction to the transmitted wave. 

The transmission coefficient 2 T and the reflection coefficient R are the squares of the limiting 
ratios of the wave amplitudes. Thus 

r=e- 2 «|l + ei(c 3 )|- 2 , (2.11) 

/?=|l + e 1 (c 2 )| 2 |l + e 1 (c 3 )|- 2 , (2.12) 

where ([1], (5.17)) 

le^c,)! *£ h exp {t CuCj (F)}-i (, = 2, 3), (2.13) 

the variations being evaluated along f -progressive curves. 

Alternative, and sometimes more effective, formulas for T and, especially, R, can be derived 
with the aid of the energy-conservation eq ([8], p. 1066; [9], sec. 4.4) 

T+R=l. (2.14) 



2 Or "tunnelling probability" [7]. 
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(This relation may be verified analytically by means of Theorem 4 of [10] and the Wronskian prop- 
erty of the solutions of the differential equation.) Substitution of (2.11) and (2.12) in (2.14) yields 



1 1 + ci(c 3 ) | 2 — 1 1 + €!(c 2 )| 2 = e" 2 «. 



He 



|l + ei(c 2 )| 2 + e- 2f * 



R = l- 



li + ^te)! 2 



l + eifoOl 2 l + e!(c 2 ) 2 + e- 



(2.15) 



(2.16) 



(2.17) 



Remark. It is interesting to observe that we have not used the connection formula for the 
turning point at z = -\- a. In contrast, if the problem of this section were attacked by real-variable 
theory ([1], sec. 9) connection formulas at both — a and + a would have to be used, and the resulting 
approximate solutions of the differential equation then matched in the interval — a < z < a; see, 
for example, [5, 6]. 3 The error analysis obviously would be lengthier and lead to more compli- 
cated bounds. 

3. Example 4 

Consider the modified Weber equation 

d 2 wldz 2 = (a 2 — z 2 )w (a > 0). (3.01) 

We take f(z) = a 2 ~z 2 and h(z) = 0. Then 

f(z) = \ {a 2 -t 2 yl 2 dt=hTTa 2 -ha 2 cos" 1 (z/a)+ iz(a 2 -z 2 y' 2 . (3.02) 

J -a 

The principal subdomains in the z-plane are shown in figure 2 (compare [12], sec. 3). The 

3 A real-variable method which avoids this matching has been suggested by Miller and Good [4] and by Pike [7], but it has yet to be placed on a firm mathematical 
foundation. 

4 For an account of some direct numerical calculations of transmission coefficients to high precision see [11]. 
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FIGURE 2. z-plane: principal subdomains. 



Figure 3. gi-plane. 
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corresponding map of the upper half of the z-plane on the £ r plane is indicated in figure 3; the map 
of the lower half is of course the image of figure 3 in the real axis. 

The reference points d, c 2 , and r :i are taken at infinity on outward-drawn rays from the origin 
in the first, second, and third quadrants of thez-plane, respectively. With the aid of the £ r map, 
we see immediately that hypotheses (i) and (ii) of section 2 are satisfied. Theorem 4 of [1] shows 
that hypothesis (hi) is satisfied, and also that €i (c 2 ) = 0. From the definition (2.06) we have a=\ira 2 . 
Hence from eqs (2.16) and (2.17), we derive 



1 + e- 



R = 



1 



1 + e" 



(3.03) 



These results may be verified analytically with the aid of asymptotic expansions and connec- 
tion formulas for the Weber functions given, for example, in [13]. The actual expression for W\{z) 
is found to be, in Miller's notation, 

wi(z) = 2-( ia2+ »l*a ia2 l 2 exp (-| rra 2 -^ ia 2 -^ Up\ Ur^W (| a 2 , V2z\ + ik^W f| a 2 , - V2z)j. 

(3.04) 

It is noteworthy that our theory 5 proves sufficiently powerful to yield exact expressions for T and 
R. This example, however, is unique in this respect. 

4. The Approximate Harmonic Oscillator 

The conditions assumed in this section are those stated in the opening paragraph of section 2, 
except that f(z) is negative when — a < z < a. 



We distinguish the auxiliary functions and principal subdomains associated with the turning 
point at z = —a from those associated with z=+a by the addition of tildes. Thus 



J (I J -(I 



(4.01) 



Three principal curves emerge from —a and three from -ha, one of which in each case is the join 
of —a and -ha; see figure 4. Let Di, Di be the principal subdomains which include the points 



5 This is also true of the theory of [14], section 9.1. 




c 3 = c 2 



Figure 4. Typical principal subdomains for the approxi- 
mate harmonic oscillator. 
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Cl = + oo and C\ = — oo ? respectively. Then clearly D 2 = D 3 and D 3 = D 2 . We take c 2 , c 3 to be 
points at infinity on conjugate ^-progressive curves lying within D 2 , D 3 , respectively, and ? 2 = c 3 , 
< r, 3 == c 2 . We assume that ^-progressive curves can be found in D linking c 2 and c 3 with -f °° and 
— oo, and also with each other, and that |Re £(z)| — > °° as z tends to any of the Cj or q. 
The fundamental connection formulas for the two turning points are given by 



{1 + € 2 (c 3 )}w 1 (z) + e~ 2 ^ 3 {l +€i(ci)Wz) + e 2 ^ 3 {l + €i(c 2 )}w 3 (z) - 0, 



(4.02) 



{ 1 + € 2 (c 3 )}wi(z) + e- 2 ^ 3 {l + € 3 (ci)}^ 2 (z) + e^{l + €i(c 2 )}w 3 (z) = 0. 



(4.03) 



Clearly w 3 (z) is a multiple of w 2 (z), and w 2 (z) a multiple of w 3 (z). To fix these relationships we suppose 
that gi (z) andgi(z) are the branches of the function {/(z)}" 1/4 which take their principal values when 
z > a and z < — a, respectively. We then find that if z E D 2 



6(z) = ia +| 3 (z), * 2 (z) = c- w ^ 3 W , 
where a is again defined by (2.06). Therefore 

ii z (z) = e (7ri/6)+ia w 2 (z) , e 3 (z) = e 2 (z). 



(4.04) 



(4.05) 



Substituting these and the corresponding results for w 2 (z) and e 2 (z) in (4.03), we obtain the connection 
formula for w x (z), w 2 (z) and w 3 (z). Elimination of w 3 (z) by means of (4.02) then yields 



Wi(z) : 



l + €i(c 3 ) jl±jl(^)3il±jl(£^ e 2ia 

l + e 2 (c 3 ) {l + € 2 (c 3 )}{l + €i(c 2 )} 



.„ , , 1 + €i(c 2 ) • 



ie ia wi{z). 



(4.06) 



The eigencondition or "quantization condition" is that Wi(z) be a multiple of ^i(z). Accord- 
ingly, this is expressed by 



^_ {l + 6i(c 3 )}{l+6i(c 2 )} 
{l + €i(c*)}{l + ii(c 8 )} 



(4.07) 



13 5 
Neglecting the error terms, we obtain at once the well-known approximations & = — 77, - 7r, — 7T, . 

To assess the errors in these approximations, we write 

ei(c 3 ) = pe^, ei(c 2 ) = pe-*, e fa) = Pe % ifa) = pe~ r % 

where p, p > and q, q are real. Then 

In 



{l + ei (C3)}{l+il(C2)} 
L{l + €i(c«)}{l+ii(Cii)}J 



= 2i(co + a)), 



where co, d> are real, and 

|o>| =^ sin -1 p ^ 2 7rp, |w| ^ sin -1 p ^ 277-p, 

provided that p, p^ 1. Hence the eigencondition is given by 

a=(M+^-h^)7r (7i = 0, 1, 2, . . .), 
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(4.08) 



(4.09) 



(4.10) 



(4.11) 



where /8 = (to + w)/ir and is real. From (4.08), (4.10) and [1], (5.17), we obtain the desired bound 

|0| « \ exp {f »,eJF)}+\ exp {*-.,*&)} ~\, (4.12) 

where F is any continuous branch of the integral defined by (2.02) above, and the variations are 
taken along ^-progressive paths. 

When the eigencondition is satisfied, eq (4.06) reduces to 

w i (z) = (-)"Y^ »«ito. (4.13) 

where 

|8|* J exp {t x , r2 (F)}-i, |fi| ^ J exp {t _ x , Cl (F)} - i. (4.14) 

Thus on the real axis the eigensolution is given by 

w 1 (z) = {l + e 1 (z)}f- 1 l 4 exp (- jV 2 ^), (4.15) 

where 

|€i(z)| ^exp {t X ,,(F)}-1 (z>a), (4.16) 

and 

M ; 1 (2) = (-)»(l+8)(l + 8)-'{l+iiW}/- ,/4 exp (-('"fl^dt), (4.17) 

where 

|€i(2)| ^exp {f _ X ,*(F)}-1 (z<-a). (4.18) 

The fractional powers in (4.15) and (4.17) are positive, and the variations in (4.16) and (4.18) are 
taken along the real axis. 

Remark. Let f(z) and h(z) be even functions of z, and assume that C2, o,\ can be taken as the 
points ±£°o. Then the bound (4.12) and eq (4.13) reduce to 

|0|^4 exp{f x , /x (F)}-i, (4.19) 

wi(z) = H n u\(z) = {-) n w x (- z). (4.20) 

5. Examples 



and 



(i) Weber equation. 



d 2 w/dz 2 = (z 2 - a 2 )w (a > 0). (5.01) 

We take f(z) = z 2 ~a 2 and h(z) = 0. Then 

f(z)= r U 2 -a 2 j 1 / 2 ^ = ~2zU 2 -a 2 ) 1 / 2 -^a 2 In {z + (z 2 - a 2 ) 1 / 2 } + la 2 In a. (5.02) 

287 




Figure 5. Principal subdomains for the Weber equation. 



Figure 5 illustrates the principal subdomains. We take C2, c 3 to be the points ±i°o. From [1], 
theorem 4, we see that €1(02) — €1(02) = 0. Hence in (4.08) p=p~0. Therefore f3 vanishes and 
a = {n+ 2)77, exactly.^ Using (2.06), we see that this implies a 2 = 2n-\-\. The eigensolution of 
course obeys (4.20). These results accord with the known analytical properties of the Weber 
functions [13]. 

(ii) Large parameter. Let/(z) = u 2 fo(z), where u is a large positive parameter, and the functions 
fo(z) and h(z) are independent of u (though this restriction could be eased without significantly af- 
fecting the following analysis). Then 



a = ua , F(z) = F (z)/u, 



whe 



OLq 



= J_* l/ol 1 ' 2 ^, Fo(z)= j {j^-jyi ;S (tt/< 



fl 14 dz 1 \f\ 



and are independent of u. Write 



dz, 



1 



X (a)=-exp {a- 1 ^ x , r2 (F )}+-exp {u~^ '.., C2 (F )} 



(5.03) 



(5.04) 



(5.05) 



Then the eigencondition (4.11) is expressed by the following equation for a: 

l//(w) = UCto — (n -h 2 ) 7T — 7Tj8(^) = 0, 

in which 

|/3(«)| « x(«)- 



(5.06) 



(5.07) 



For large «, x(w) =0(u '). Therefore for large values of the positive integer n 

u=oq\n+\ )TT + 0(n- 1 ). 



(5.08) 



6 Compare also [14], page 109. 
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In order to determine precise bounds for the O term in the last equation, consider the function 
X(&). This tends to infinity as a^ 0, and to zero as u— » «>, and is monotonic strictly decreasing 
when u is positive. Let u be the root of x( u o)= 2 . Then provided that nir/a Q > u , we have from 
(5.06) and (5.07) 

i//(/i7r/a ) = - i 77 - 7Tj8(/i7r/ao) < 0, (5.09) 

and 

i/f{(/i + l)77/a } = %7T-TT(3{(n+ 1)tt/o:o} > 0. (5.10) 

Therefore at least one eigenvalue exists in the interval nrr/ao < u < {n -h l)7r/ao. To delimit it in 
a shorter interval, let 

a = (n+ i+i;)ir/ao, (5.11) 

where |v|< 2. Then substituting in (5.06), we obtain 

1; = £{(»+ J +t;)7r/oo}. (5.12) 

Hence \v\ < xi^^l^o), that is, 

H <|exp l^-r x>C2 (Fo)} + |exp f— ^-*,c 2 (Fo)]-|. (5.13) 

4 [n77 J 4 IniT > J z 

Relations (5.11) and (5.13) constitute the precise form of the eigencondition (5.08), and are 
valid when n > OqUqItt. The eigensolution satisfies (4.13), where 



|S| *£- exp • 



~ ^x,c 2 (Fo)]-|, |6|^5 exp {^ •/ _ x , C2 (F o) ]-l. (5.14) 

M77 J Z ^ [7177 J Z 



The above proof shows that there is at least one eigenvalue (5.11) which satisfies (5.1.3). In 
order to establish that there is exactly one eigenvalue of this form, we need to investigate the sign 
of the derivative /3'(u) in the interval considered. 7 This takes us back to fundamentals. We have 
to consider the Liouville-Green approximate solutions of the differential equation 

d 2 wldz 2 ={u 2 +f(z)}w (5.15) 

(compare [1], sec. 2) as functions of u, and seek bounds for the ^-derivatives of their error terms. 
The analysis is somewhat lengthy, and we record only the result: 
Let 

v/ , W 1 + k, le 1 1+k 

X(u) = - - + - r oo, (5.16) 

d — e' a o~-e l u 

where 

l^w^^eJiFoX l = u- l f _ X ,,. 2 (F () ), (5.17) 

K = (2e cos y)~\ k= (2e cos 7)" 1 , (5.18) 



7 An alternative, but apparently no easier, approach would be to consider complex values of u and apply Rouche's theorem. 
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and y, y are the maximum angles of slope with the real £-axis of the f-maps of the progressive paths 
linking o° to c 2 , and — °° to c 2 , respectively (|y|, | y| < 2 if). Then exactly one eigenvalue satisfies 
(5.11) and (5.13), provided that n ^ a u^7T where U\ is the largest positive zero of X{u). 

Some simplification can be effected by using upper bounds for the zeros of x( u )~ 2 andX(u). 
In this way we may verify that both of the conditions n > (XoUqItt and n ^ a Uilrr are certainly ful- 
filled if n exceeds each of the three quantities 

it In 1 77 In 2 77 ' 
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